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1 . 0 Introduction 


During the past feu years there has bean an increasing interest 
in the effectiveness and efficiency of ice protection systems. To 
address these concerns, a grant from NASA Lewis Research Center 
(NASA LeRC) was awarded to the University of Kansas in May of 1980 
entitled "An Assessment of General Aviation Deicing Systems on Con* 
temporary Airfoils.” Two ice protection systems were evaluated on 
two general aviation wing sections over a period of two years as 
part of this grant. The inflatable tube pneumatic boot was one of 
these systems. The second system, originated in England in the 
1940s, was the fluid ice protection system. This system incorporates 
a porous leading edge panel through which a glycol-water solution 
exudes on to the wing to lower the freezing temperature of the 
resulting mixture with the impacting water. The results of the 
fluid system evaluation are presented in ref. 1. 

As a result of these tests, the fluid ice protection system has 
rapidly increased in popularity with U.S. aircraft manufacturers. 
Today, for example, the Cessna Citation IISP is flying with the 
system as factory installed equipment. Beechcraft's Starship I and 
the Old Man's Aircraft Company OMAC-I will also be offering the 
system as standard equipment, flight test and certification are now 
underway for a Supplemental Type Certificate (STC) on the Cessna 
206. In Great Britain, the Beech Duchess, Piper Seminole, and the 
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Grumman American Cougar have already been approved for flight in 
known icing conditions with a fluid system. Tests are now underway 
in Great Britain for approval of the fluid system on the Beech Baron 
and the Partenavia P-68B. 

With the increasing popularity of the fluid ice protection 
system, it has become necessary to establish a broad data base of 
®®ti ice flow rates for a variety of icing conditions, wing sec~ 
tions, and types of porous panels. Also, there is a need for a 
verified method of predicting anti-icing flow rates. Such a method 
would reduce the necessity of expensive wind tunnel tests and would 
be a valuable tool for the airplane designer. 

The research presented in this report was supported in full by 
a follow-on grant from NASA LeRC to the University fo Kansas en- 
titled "A Joint KU/NASA Graduate Research Study." The objectives of 
this research effort are as follows: first, to increase the ex- 
perimental data base of minimum required anti-icing flow rates 
throughout a broad range of icing conditions on two different 
general aviation airfoils; second, to test the validity and 
reliability of two methods of predicting anti-ice flow rates; 
third, to test a new laser drilled titanium porous panel; and 
finally, to evaluate a composite porous panel. 

This report is presented in four sections. The first section 
presents the test facility, test conditions, and the test procedures 
used. The second section is a brief discussion of the aerodynamic 
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effects associated with ice accretions. The third section presents 
the anti-ice flow rates obtained through wind tunnel tests in the 
NASA LeRC Icing Research Wind Tunnel. The final section is a 
presentation of the prediction methods, and a comparison between the 
predicted and actual flow rates. 

2.0 Description of Test Facilities and Equipment 

2. 1 Icing Research Tunnel 

The NASA Lewis Icing Research Tunnel (IRT) is a closed cycle, 
refrigerated wind tunnel with a rectangular test section 1.33m (6 
ft) high by 2.74 m (9 ft) wide by 6.1 m (20 ft) long (Fig. 2.1). 
Maximum tunnel airspeed is approximately 134 m/s (300 mph) . A 
natural icing cloud is simulated by injecting a water spray upstream 
of the test section from a multi-nozzle spray bar system. 

The area of interest on a test model is confined to the region 
in the center of the test section, where the icing cloud is most 
uniform, covering a cross-sectional area of approximately 0.9 m (3 
ft) high by 1.5 m (5 ft) wide. The liquid water content (LWC) of 
the cloud can be varied from about 0.5 to 2.4 g/m 3 with volume 
median droplet diameters (d) in the range from 10 to 20 microns. 
The tunnel total air temperature can be regulated from ambient to 
below -28.9 degrees C (-20 degrees F). 
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From outside the tunnel test section, steam and electrically 
heated windows allow visual observation and documentation of the 
test while in progress. Between icing sprays, tunnel operators and 
engineering test personnel, appropriately dressed with cold weather 
gear, are able to enter the tunnel test section. This access to the 
test model permits manual cleaning of ice and frost off the wing 
model, documentation of ice shapes, close-up photography, and 
modifications to the test model. 

A translating wake-survey probe was used to measure the section 
drag coefficient, c^, of the test ncjdel. The probe, which had a 
stagnation and static pressure orifice, could be retracted behind a 
wind screen. When the airfoil was exposed to the tunnel icing 
cloud, the probe was retracted. After the icing cloud was turned 
off, the probe was inserted into the airstream and the translating 
drag survey was made. The probe was located about one chord length 
downstream of the airfoil at midspan. 
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Figure 2.1. - Schematic Planform of Icing Research Tunnel 



2.2 Wing Models 


This report contains results and analysis of tests conducted on 
two different general aviation wing models, designated as Wing Model 
A and Wing Model B. These wing models, having different airfoil 
sections, were tested independently. Results and conclusions ob- 
tained on one wing model cannot be generalized to other wings. 

Wing Model A was an inboard section of a full scale wing taken 
from a single engine general aviation aircraft. The wing model had 
a constant chord length of 1.63 m (64.0 in) and a modified leading 
edge NACA 2412 airfoil as depicted in Fig. 2.2. Table I presents 
the coordinates for the airfoil. The model was not twisted or 
swept. The wing model was fastened to the tunnel floor turntable 
through brackets bolted to the primary internal wing structure. It 
was supported at the ceiling through a pivot. A clearance of 1/2 
inch was allowed between the model and the ceiling of the 6 foot 
high test section. 

Wing Model B also was a full scale wing taken from a single 
engine general aviation airplane. The original wing tapered from a 
NACA 64^A215 airfoil at the root (WS 0) to a NACA 64^A412 airfoil at 
the tip (WS 216). The wing incorporated a leading edge modification 
designed by Raymond Hicks (ref. 2,3). This modification increases 
the maximum lift coefficient and the lift/drag ratio. The thickness 
over the forward 30 percent of the upper surface and the leading 


I 


edge radius is increased as shown in Figure 2.3 on a NACA 64^A412 
airfoil . 

The wing section was fastened to the tunnel floor turntable 
using the spar fittings that are designed to attach the wing to the 
fuselage of the airplane. The centerline of the tunnel was at WS 58 
of the original wing. Table II gives the airfoil coordinates at WS 
58, where the chord is 1.61 m (63.25 in). The chord tapered 1.1 
inches per foot of span, and the wing was twisted 0.167 degrees per 
foot of span (washout). Figure 2.4 shows Wing Model B in the IRT 
test section. 
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MODIFIED NACA 2412 



Figure 2.2. - Airfoil Section of Wing Model A 


NACA 64^412 

• *' ** PropoMd Modification 
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Figure 2.3. - Hicks Modification on a NACA 64.A412 Airfoil 
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Airfoil Coordinates of Wing Model A in 
Percent of Chord. 


Upper 

Surface 

%x/c 

XY/C 



0.00 

0.00 

1.25 

2.03 

2.50 

2.93 

5.00 

4.06 

7.50 

4.97 

10.00 

5.64 

15.00 

6.55 

20.00 

7.36 

25.00 

7.71 

30.00 

7.72 

40.00 

7.90 

50 . 00 

7.38 

60.00 

6.41 

70.00 

5.24 

80.00 

4.02 

90.00 

2.71 

95.00 

2.08 

100.00 

1.35 


Lower Surface 


%x/c 

XY/C 

0.00 

0.00 

1.25 

-3.44 

2.50 

-4.06 

5.00 

-4.06 

7.50 

-4.02 

10.00 

-4.00 

15.00 

-4.06 

20.00 

-4.29 

25.00 

4.29 

30.00 

1 

t— * 
00 

40.00 

-3.84 

50.00 

-3.39 

60.00 

-2.91 

70.00 

-2.23 

80.00 

-1.17 

90.00 

.09 

95.00 

.72 

100.00 

1.35 
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Table II : Airfoil Coordinates of Wing Model B at Centerline 

of Tunnel (WS 58) in Percent of Chord. 


Upper Surface 


Vi/c 

Vi/C 




.000 

- .704 

.015 

- .250 

.648 

.791 

1.138 

2.372 

2.055 

3.477 

3.953 

4.941 

6.324 

6.008 

9.486 

6.735 

11.352 

7.036 

13.439 

7.502 

22.024 

7.565 

24.996 

7.581 

30.126 

7.597 

34.783 

7.543 

39.428 

7.426 

44 . 409 

7.110 

49 . 387 

6.591 

54.360 

5.891 

59.331 

5.047 

62.111 

4.526 

99.744 

-2.606 


Lower Surface 
Vi/C XY/C 


.000 - .704 

.335 -1.474 

•723 -1.858 

1.216 -2.193 

2.451 -2.760 

4.926 -3.545 

7.407 -4.130 

14.223 -5.371 

19.197 -5.395 

24.175 -6.359 

29.157 -6.658 

34.142 -6.816 

39.129 -6.870 

44.122 -6.718 

49.115 -6.449 

54.111 -6.114 

99.741 -2.794 
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2.3 Ice Protection System* 

2.3.1 System Description 

The fluid ice protection system tested consists of a porous 
panel attached to the leading edge of a wing, and a pumping system 
that distributes a glycol based freezing point depressant fluid from 
a tank to the panel. The rate of fluid exuding from the porous 
panel onto the surface of the wing and the severity of the icing en- 
counter determines the type of ice protection possible. The system 
is capable of anti-icing by creating a mixture of the glycol fluid 
with the impacting water with a freezing temperature below the local 
temperature to prevent the formation of ice on the leading edge. 
The system is also capable of deicing by dissolving the bond between 
the accreted ice and the panel skin until aerodynamic forces carry 
the ice accretion downstream. 

Most ice protection systems operate in either an anti-ice, or a 
deice mode. However, the fluid ice protection system is capable of 
operating in, what is often referred to as the natural-deice mode. 
In this mode the flow of glycol is lower than the anti- ice flow 
rate, but still high enough that no permanent ice accretions form on 
the wing leading edge. In the natural-deice cycle, temporary ice 
accretions ranging from small ice patches to narrow bars form near 
the stagnation point, before being carried downstream after several 
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seconds to a few minutes. This build-shed cycle then repeats, 
without forming any permanent accretions of ice. 

A benefit of the fluid ice protection system is the runback of 
the glycol solution onto t u e upper and lower wing surfaces behind 
the leading edge panel. This runback is often sufficient to pre- 
vent, or remove, the formation of any aft frost which might exist on 
the wing surfaces. 

2.3.2 Porous Leading Edge Panels 

Three types of porous leading edge panels were tested; stain- 
less steel mesh, laser drilled titanium, and composite. The stain- 
less steel type panels were tested on both wing models. The 
titanium panel was tested only on Wing Model A. The composite panel 
was tested only on Wing Model B. The stainless steel and composite 
panels were riveted onto the original wing leading edge skin. The 
titanium panel was attached with countersunk screws. Each of the 
panels tested were divided into three independent sections into 
which the flow rate of glycol could be controlled independently. 
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2 . 3 . 2 . 1 Stainless Steel Mesh Panel 


The stainless steel mesh panel consisted of two layers of wire 
cloth that were rolled, sintered, and then finish- rolled. The wire 
cloth was manufactured from an 18-8 austenitic stainless steel and 
nominally has 110 by 20 wires per inch. The two layers of wire 
cloth are laid 90 degrees with respect to each other. Figures 2.5 
a-c show a front view, back view, and a close-up view of the stain- 
less steel mesh panel. 

Inside the fluid reservoir, behind the stainless steel mesh 
skin, is a thin polyvinylchloride sheet. The purpose of this 
material, whose porosity is much lower than that of the stainless 
steel, is to increase the pressure in the resevoir sufficiently to 
distribute the glycol evenly over the entire active portion of the 
panel, regardless of the chordwise pressure distribution, which 
changes with angle of attack. 

The maximum thickness of the leading edge panel was 3.2 mm. 
The panel contoured to Wing Model A section had an active porous 
width of 6.9 cm. The length of the inboard and outboard sections 
was 29.2 cm long, and the middle section was 40.6 cm long. The 
panel installed on Wing Model B had an active width of 8.7 cm, with 
the outboard and inboard sections 29.2 cm long, and the middle sec- 
tion 40.6 cm long. 
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2. 3. 2. 2 Drilled Titanium Panel 


The titanium panel is similar to the stainless steel panel ex- 
cept for a laser drilled titanium outer skin perforated with 
microscopic holes through which the glycol fluid passes. The holes 
are nominally .0025 inches in diameter, with .035 inch separation 
between centers. This yields approximately 800 holes per square 
inch. The panel had an active width of 7.6 cm. The inboard and 
outboard sections were 29.9 cm long, and the middle section was 41.3 
cm long. Figures 6 a and b show a back view and a side view of the 
drilled titanium panel tested on Wing Model A. 

2.3.2. 3 Composite Panel 


The porous composite panel, tested only on Wing Model B, had an 
outer skin made with a resin starved matrix laminate with successive 
layers of Kevlar49 - 3D woven Eglass - Kevlar49. Inside the fluid 
cavity were layers of polyvinyl fluoride, which is intended to 
distribute the flow of glycol uniformly chordwise regardless of the 
pressure distribution changes that occur as angle of attack changes. 
The original panel tested in 1981 (ref. 4) contained only 1 layer of 
the polyvinyl fluoride material. The second composite panel, 
discussed in this report, contained an additional layer of this 
filter material. The backing plate was constructed of an ordinary 
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boat glass epoxy layup. The width of the porous region was 8.7 cm. 
The inboard and outboard section were 29.2 cm long and the middle 
section was 40.6 cm long. The maximum thickness of the panel was 
3.2 ran. Figures 2.7 a and b show the front view and back view of 
one of the three sections of the composite panel. The front view 
illustrates the distribution of fluid exuding out of the panel 


during a static test. 
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(a) Front View 



(b) Back View 

Figure 2.7. - Composite Porous Leading Edge Panel 


2.3.3 Freezing Point Depressant Solutions 


Two glycol based fluids were used during the tests presented in 
this report. The first fluid's composition (TKS80) , as determined 
by the NASA Lewis Chemical Laboratory, was 18.8% de-ionized water 
and 81.2% mono-ethylene glycol. The second fluid's composition 
(AL5) was determined to be 12.2% de-ionized water, 3.1% iso- 
proponal, and 84.7% mono-ethylene glycol. Both fluids have specific 
gravities of 1.09 g/cm 3 . 

The tests reported in reference 1 were conducted using only the 
TKS80 fluid. All the tests presented in this report used the ALS 
fluid. In addition, the TKS80 fluid was used during the stainless 
steel panel test on Wing Model A to determine any differences in 
performance between the two fluids. Figure 2.8 presents the 
freezing temperatures (T^) of the solutions as they are diluted with 
water. The temperatures are plotted against the mass fraction of 
glycol plus proponal , G, in the mixture. At mass fractions between 
approximately 65 and 80 percent the mixture turns to a gel instead 
of a solid as it is chilled. 
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2.3.4 Fluid Supply Pump 


The flow rate to each of the three sections of a porous leading 
edge panel was controlled and metered independently by a flight- 
worthy pumping system with three variable positive displacement 
pumps. The overall pumping package as shown in the schematic of 
Figure 2.9 consisted of a transfer pump, fluid filter, air filter, 
fluid tank, control panel (not shown), and the main metering pump. 
Figure 2.10 shows the actual pump system apparatus used during the 
wind tunnel tests. 

The main fluid tank is filled using a transfer pump which pumps 
the fluid from the original container through a filter and into the 
tank. The tank outlet is connected to a small air filter to prevent 
the ingress of dirt. A float type contents transmitter is installed 
in the tank. The control panel provides control of and indicates 
all pump functions along with an indication of the level of fluid in 
the tank. 

The main metering pump consists of three separate pumping 
cells, driven by a 28 volt DC motor through a toothed belt reduction 
drive train. Each cell consisted of a piston actuated by an eccen- 
tric cam. A stroke limiting device, connected to a digital position 
readout is used to provide an adjustable flow rate. Pressure relief 
valves, bolted on the outlet of each cell, are set to open at ap- 
proximately 95 lbs/sq. inch to protect the porous panel from exces- 
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sive pressures. A solenoid valve connected to each pump cell 
provides a means of selectively turning each outlet on or off as re- 
quired . 

The fluid ice protection system, including the stainless steel 
and titanium panels, TKS80 fluid, and the pumping system were 
developed and manufactured by T.K.S. (Aircraft De-Icing) Limited, of 
Great Britain. The North American representative of T.K.S. is Kohl- 
man Aviation Corporation in Lawrence, Kansas. The composite porous 
panel was developed and manufactured by Fiber Materials Incorporated 
in Biddeford, Maine. The AL5 fluid, marketed in the United States 

and approved for use by T.K.S., was supplied by Canyon Industries of 
Phoenix, Arizona. 
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Figure 2.9. - Schematic of Fluid Supply Pump System 
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(a) Front View 

Figure 2.10. - Fluid Pump System Apparatus Used During 
Tunnel Tests. 
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3.0 Test Techniques 


3.1 Tunnel Conditions 


Test airspeed conditions of 49.2 m/s (96 knots), 74.6 m/s (145 
knots), and 90.2 m/s (175 knots) were selected based on the best 
climb and cruise speeds of aircraft on which the airfoil sections 
tested are used. Operating envelopes of liquid water content and 
drop size, dependent upon tunnel airspeed, were plotted for the air- 
speeds chosen. From these tunnel operating envelopes several com- 
binations of LWC and drop size were chosen for testing, including 
the least and worst icing severity conditions of each envelope, as 
plotted on Figures 3.1 a .. The type of ice (i.e. glaze, rime) that 
formed on the wing model depended primarily on the tunnel total air 
temperature (T) . To accrete primarily glaze ice, the air tempera- 
ture was set at -3.9 degrees C (25 degrees F) ; and to accrete 
primarily rime ice, it was set at -15.0 degrees C (5 degrees F), or 
at -23.3 degrees C (-10 degrees F). Figures 3.2 and 3.3 present the 
continuous maximum and intermittent maximum icing condition curves 
specified in Federal Aviation Regulations, Part 25 (ref. 5). 
Figures 3.4 and 3.5 overlay the IRT operating envelopes at 96, 145, 
and 175 knots with the continuous and intermittent maximum en- 
velopes . 
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Figure 3.1. - IRT Operating Envelopes and Test Points 
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Figure 3.2. - Federal Aviation Regulations Part 25 (ref. 5 ) 
Continuous Maximum Icing Conditions. 
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Figure 3.3. - Federal Aviation Regulations Part. 25 (ref. 5 ) 
Intermittent Maximum Icing Conditions. 
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3.2 Test Procedures and Analysis 


3.2.1 Aerodynamic 

Aerodynamic data were recorded for Wing Model A in a clean con- 
dition (no ice or frost) and with ice on the wing. The clean wing 
test consisted of determining the location of stagnation points 
about the leading edge of the airfoil throughout a broad range of 
angles of attack. This was accomplished with the porous panel by 
placing 18 static ports on strip-a-tube laid chordwise at 1/2 inch 
intervals. The strip-a-tube was attached to the wing surface with 
adhesive tape. After stabilizing at the appropriate tunnel air- 
speed, the peak surface static pressures were recorded at each angle 
of attack. 

The second type of aerodynamic test was performed with the ice 
protection system turned off, allowing ice to accrete. The 
procedure for this test was as follows. 

a) start with an ice and frost free wing 

b) set velocity, angle of attack, and temperature 

c) set icing cloud conditions, LWC and d 

d) turn on icing cloud 

e) after desired icing encounter, turn cloud off 

f) take drag wake survey 
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g) wipe off any surface frost behind primary ice formation 

h) take a second drag wake survey 

i) clean off ice in preparation for next run 

3.2.2 Anti-Ice 

The procedure for obtaining minimum required anti-ice fLuid 
flow rates was as follows. 

a) set velocity, angle of attack, and temperature 

b) set icing cloud conditions, LWC and d 

c) turn on glycol flow rate to center panel at 

a rate in excess of that required for anti-icing 

d) turn on icing spray 

e) slowly reduce glycol flow rate rate until 
minimum anti-icing is achieved 

f) turn icing spray off and prepare for next run 

Initially at a given flight condition, the glycol fluid flow 
rate was set to be well above that required for anti-icing. The 
flow rate was then reduced in steps, allowing about 20 seconds for 
the system to stabilize, until small flecks of ice began to appear 
on the leading edge in the vicinity of the stagnation point. At the 
anti-ice flow rate small ice flecks would just begin to form before 
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being swept downstream in only a few seconds. A glycol flow rate 
slightly lower would cause the ice flecks to persist, gradually 
growing into larger patches before being shed from the wing. Often 
the first sign of ice would be a small fleck, on a rivet head in the 
porous region of the panel near the stagnation location. 

The center sections of the porous panels, in the middle of the 
tunnel icing cloud, were used to obtain the minimum anti-icing fluid 
flow rates. However, several aspects made it difficult to obtain 
these flow rates with a high degree of accuracy. 

First, at a given glycol flow rate, regions of the center panel 
would be anti-icing, while at a different spanwise location, the 
panel would be operating closer to the natural-deice mode. It was 
difficult to determine the appropriate anti-ice flow rate that was 
characteristic of the entire center panel. Second, on many of the 
icing conditions, the transition from anti-icing to deicing was 
subtle, as it occurred over a broad range of glycol flow rates. It 
was therefore difficult to determine the unique minimum anti-ice 
flow rate. Finally, it was difficult to observe the small ice par- 
ticles forming and shedding on the wing from a distance of five 
feet, while looking through an icing cloud. This problem was par- 
tially alleviated by sighting through a short range telescope with a 
field of view of approximately 4 inches in diameter on the leading 
edge. 
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Because of the subjectiveness involved in determining the anti- 
ice flow rates, each unique icing condition was tested from one to 


five times. 

An average flow rate at 

each icing condition was 

then 

calculated . 

It is 

believed that 

this average 

value is 

more 

representat 

ive of the 

true anti,-ice 

flow rate, than 

a single 

test 


point. Because of the subjectiveness involved in obtaining these 
flow rates, a second observer, a NASA engineer in the Aircraft Icing 
Research Section, independently judged flow rates during the test 
series on Wing Model B with the stainless steel panel. Equal weight 
was given to the two observations in determining the average flow 
rate. 

3.2.3 Deice 

The test procedure for determining the time required to shed an 
ice accretion at various icing conditions and glycol flow rates was 
as follows: 

a) set velocity, angle of attack, and temperature 

b) set glycol flow to inboard and outboard panels 
sufficient for anti-icing 

c) set icing conditions, LWC and d 

d) turn spray on to accrete ice on center panel 

e) turn spray off after desired icing encounter 
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f) turn glycol flow on to center panel 

g) record elapsed time after beginning glycol flow 
to center panel until ice cap sheds 

h) clean ice and frost off wing in preparation for 
next run. 

4.0 Test Results and Discussions 

4. 1 Aerodynamic 

Figures 4.1.1 and 4.1.2 present the stagnation point locations 
for angles of attack ranging from -2 to 12 degrees on Wing Model A 
and B, respectively. Stagnation locations for angles ranging from 0 
to 11 degrees fall within the porous region of the leading edge 
panel installed on Wing Model A. On Wing Model B, stagnation loca- 
tions for angles from approximately -1.5 to over 12 degrees fall 
within the porous region. 

As reported in reference 1, the installation of the stainless 
steel porous panel on Wing Model B adds less than .001 to the sec- 
tion drag of the wing, within the uncertainty band of the drag 
measuring system. 

In order to obtain an understanding of the aerodynamic penal- 
ties associated with ice accretions on a wing leading edge, a brief 
test was conducted on Wing Model A with the ice protection system 
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inoperative. The test runs were conducted at a velocity of 145 
knots and an angle of attack of 2.1 degrees - 

Glaze ice, accreted at 25 degrees F, is characterized by the 
typical double horn ice formation. As a result, the drag penalties 
associated with an extended glaze icing encounter can be severe. 
Figures 4.1.3 a and b illustrate the percent increase in drag coef- 
ficient above the clean wing drag, c^, for glnze icing at two li- 
quid water contents. Figures 4.1.4 and 4.1.5 present pictures of 
the corresponding ice shapes . 

Rime ice, accreted at 5 degrees F, conforms more to the 
original airfoil contour as the impinging water droplets freeze on 
impact The resulting drag increases associated with rime icing are 
presented in Figures 4.1.6 a and b. The pictures of the cor- 
responding ice accretions are given m Figures 4.1.7 and 4.1.8. 

3 

After a 15 minute glaze icing encounter at a LWC of 1.16 g/m 

and drop diameter of 15 microns, the percent increase in drag trora 

the ice cap alone was 74 %. However, a rime ice accretion increased 

the drag only 22 %. This higher drag penalty associated with glaze 

3 

icing also occurs at the higher liquid water content of 1.55 g/m 
where a 15 minute glaze icing encounter increased the drag by 176 %. 
However, a rime ice accretion increased the drag only 73 % . 

Also plotted in Figures 4.1.3 and 4.1.6 are the percent in- 
creases in drag from surface frost. This frost, which often forms 
on the upper and lower surfaces of a wing section in che IRT, may 
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- **«*-«*<** i^iu(( cucounter. it is interesting to 

note the large increases in drag associated with the frost layer 
alone. After a 15 minute icing encounter, the additional percent of 
drag increase ranged from 50 to 100 percent above that from the ice 
cap alone. 



Figure 4.1.1. - Stagnation Point Locations on Wing Model A. 



Figure 4.1.2. - Stagnation Point Locations on Wing Model B. 


V/- 145 KNOTS T«25*F 
240. LWC* 155^/ d*2* 
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Figure 4.1.3. - Drag Increase as a Result of a Glaze Icing Encounter 
at V = 145 Knots, T = 25° F, and Angle of Attack = 2 
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15 microns. 



Figure 4.1.4. - (Continued.) 
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figure 4.1.5. - Glaze Ice Accreted at V = 145 Knots, T = 25° F, Angle of Attack 
LWC = 1.55 g/m 5 , and d = 20 microns. 





Figure A.1.5» - (Continued.) 
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(d) 15. Minute Icing Spray 
Figure A. 1.5. (Continued.) 
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Figure 4.1.7. - Rime Ice Accreted at V = 145 Knots, T = 5° F, Angle of Attack 
LWC = 1.16 q/m 3 , and d = 15 microns. 







(b) 10. Minute Icing Spray 
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(d) 15. Minute Icing Spray 
figure 4.1.8. - (Continued.) 


4.2 Anti-Ice 


Normal operation of the fluid ice protection system is in the 
anti-ice mode (Fig. 4.2.1); that is, the glycol flow rate is suf- 
ficient to prevent ice from forming on the leading edge of the wing. 

This is possible as long as the glycol-water mixture at the leading 
edge maintains a freezing temperature below the ambient air tempera- 
ture. The freezing temperature of the mixture increases as the 
ratio of the water catch rate to the glycol flow rate increases. A 
series of runs was conducted in the NASA LeRC IRT to determine the 
ninimum fluid flow rate at which anti-icing could be maintained as a 
function of angle of attack, liquid water content, and drop diameter 
on the following wing-panel configurations. 

>* 

- Wing Model A, Stainless Steel Mesh Panel 

- Wing Model B, Stainless Steel Mesh Panel 

i 

- Wing Model A, Drilled Titanium Panel r 

The anti-ice flow rates are presented in terms of specific 
fluid flow; grams of glycol fluid per square centimeter of active 
panel per minute. Tables A.l through A. 6 in Appendix A present the j 

i 

\ 

complete summary of the wind tunnel results on the three configura- j 

I 

tions. A significant portion of these results is presented in f 

i 

Figures 4.2.2 through 4.2.11. The symbols represent the actual 
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recorded d.t. as listed in Appendix A. The solid and dashed lines 
connect the average fluid flow values calculated at each icing con- 

dition. 
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Figure 4.2.1. - Anti-Iced Center Section of Stainless Steel 
Mesh Panel and Wing Model. 
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4.2.1 Stainless Steel Mesh Panel 


Tables A.l and A. 2 present the anti-ice flow rates determined 
on the Wing Model A with the stainless steel panel. Table A.l lists 
the runs using the AL5 fluid. Table A. 2 lists the runs using the 
TKS80 fluid. From these Tables, Figures 4.2.2, 4.2.3, and 4.2.4 
were plotted showing the anti-ice fluid flow versus ingle of attack, 
liquid water content, and volume median drop diameter respectively. 

A few general observations can be made concerning the results. 
First, the minimum anti-ice fluid flow rate increases as the stagna- 
tion location comes close to an edge of the porous region. This in- 
crease is noticed at both high and low angles of attack. At a high 
angle of attack, for example, the majority of the water droplets im- 
pact on the lower surface, while most of the glycol solution flows 
over the upper surface. As expected, the anti-ice fluid flow in- 
creases with liquid water content. Also, a noticeable increase in 
required flow rate occurs as the drop diameter is increased from 11 
to 15 microns. 

Tables A. 3 and A. 4 present the anti-ice flow rates determined 
on Wing Model B with the stainless steel panel. Table A. 3 lists the 
runs judged by Observer 1. Table A. 4 lists the runs judged by Ob- 
server 2. From these Tables, Figures 4.2.5, 4.2.6, 4.2.7, and 4.2.8 
were plotted showing the anti-ice fluid flow versus angle of attack, 
liquid water content, volume median drop diameter, and tunnel air- 
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speed . 

On this configuration, the anti-ice specific fluid flow does 
not appear to be as strongly affected by angle of attack. This is 
verified by referring back to Figure 4.1.2, showing that the porous 
region of the panel encompasses a broader range of stagnation point 


locations. 
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96 Knots, d = 20 Microns, AL5 Fluid 
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Figure 4.2.4. - Minimum Anti-Ice Fluid Flow Rates Versus Volume Median 
Drop Diameter on Wing Model A with Stainless Steel Mesh 
Panel . 
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Figure 4.2.4. - (Continued. 
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Figure A. 2. 6. - Minimum Anti-Ice Fluid Flow Rates Versus Liquid Water 

Content on Wing Model B with Stainless Steel Mesh Panel. 
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Figure 4.2.7. - Minimum Anti-Ice Fluid Flow Rates Versus Volume Median 

Drop Diameter on Wing Model B with Stainless Steel Mesh Panel. 
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4.2.2 Titanium Panel 


Tables A. 4 and A. 5 present the anti-ice flow rates determined 
with the new drilled titanium panel on Wing Model A. Table A. 4 
lists the rims determined by Observer 1. Table A. 5 lists the runs 
determined by Observer 2. From these Tables, Figures 4.2.9, 4.2.10, 
and 4.2.11 were plotted showing the anti-ice fluid flow versus angle 
of attack, liquid water content, and volume median drop diameter. 
The anti- ice fluid flow is independent of angle of attack as long as 
the stagnation location is within the center 50 to 75 percent of the 
active porous region of the leading edge panel. Figure 4.2.12 
presents a comparison of the required flow rate for anti-icing using 
the stainless steel panel with the flow rate required when using the 
titanium panel. There appears to be no significant difference in 
the required flow rate for anti-icing between the two panels. The 
points lying outside of the 25% boundary lines are an indication of 
the scatter in the experimental data. 
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Figure 4.2.9. - Minimum Anti-Ice Fluid Flow Rates Versus Angle of 
Attack on Wing Model A with Titanium Panel. 
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Figure A. 2.9. - (Continued. 
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Figure 4.2.9. (Continued. 






LIQUID WATER CONTENT > LWC > <^l nr» 4 
145 Knots, T = 25° F , d = 15 microns 
Figure 4.2.10. - (Continued.) 
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Minimum Anti-Ice Fluid Flow Rates Versus 
Volume Median Drop Diameter on wing Model 
A with Titanium Panel at V = 145 Roots t 
= 25° F, LWC = 1.16 g/m\ 


V* 96,145 KlvfcJS t ~ .5® 2. .1® 5.9° 1.1" 



Figure 4.2.12. - Comparison of Anti-Ice Flow Rates Between 

Stainless Steel Mesh Panel and Titanium Panel 
on Wing Model A. 


4.2.3 Composite 


A composite panel was also tested on Wing Model B. However, 
after only two days of evaluation, the wind tunnel tests were can- 
celled. Overall, the results obtained agreed with those from the 
original composite panel tests as reported in ref. 4. The flow 
rates required for anti-icing were two to three times higher than 
those required when using the stainless steel or titanium panel. 

These higher flow rates were the result of the net porosity of 
the panel being insufficient to distribute the flow of glycol 
uniformly in a chordwise direction. As the porosity of a 
distributing panel increases, the internal back pressure decreases. 
With low pressurization, the flow of glycol does not exude uniformly 
in a chordwise direction against an external aerodynamic pressure 
gradient. While the required glycol flow rate is the highest near 
the stagnation location, the actual glycol flow rate was the lowest 

there. 

Several other problems were also experienced with the composite 
panel. First, significant differences were noted between each of 
three independent sections. At a given setting, the resulting back 
pressure of each section was different. Second, the glycol flow 
rate varied from one region to the next within the same section. At 
one location, the flow rate was high, but at another location the 
flow rate was much lower. Third, the flow of glycol was able to 


- 85 - 


leak internally from one section to the next. Since the wing sec- 
tion was mounted vertically in the wind tunnel, a portion of the 
glycol pumped to the top panel section would slowly drain internally 
to a lower section before exuding. The combination of these 
problems caused the test series to be cancelled. 

It is believed that most of these problems are not a result of 
the external skin being a composite. Instead, it is felt that the 
problems are associated with the design, manufacturing, and quality 
control of the porosity of the panel which is primarily controlled 
by the internal filter material, polyvinyl flouride. Because of the 
significant weight savings and special applications for a composite, 
a further effort should be undertaken to develop a usable composite 
panel . 

4.3 Deice 


Deicing is an alternative operating mode for the fluid ice 
protection system. This mode would be required if ice were to form 
prior to turning the system on. Results of deicing tests on Wing 
Model A with the stainless steel mesh panel are shown in Figures 
4.3.1 a and b. These tests were conducted at a velocity of 145 
knots, angle of attack of 2.1 degrees, and after a 5 minute icing 
spray. The elapsed times required to shed the accreted ice cap are 
plotted versus the flow of glycol in terms of the percent of anti- 
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icing flow required at the respective icing condition. 

An ice accretion formed at 5 degrees F required less time to 
shed than an accretion formed at 25 degrees F. This is especially 
prominent in Figure 4.3.1a with a LWC =0.80 g/m 3 and drop size = 11 
microns. An explanation is that at 25 degrees F glaze ice is 
formed, which is accompanied by runback icing, thus a wide ice cap. 
At 5 degrees F, the ice freezes on impact with the wing, resulting 
in a narrower ice cap. The wider the ice cap, the more ice at the 
interface which must be melted before a shed is likely. 

At a glycol flow rate of 50 % of that required for anti-icing 
it was necessary to make a small change in the angle of attack to 
facilitate a shed at the elapsed times shown in Figure 4.3.1. Shed 
times of three minutes or less are possible for the icing conditions 
shown with a flow rate equal to that required for anti-icing. Only 
marginal improvement in shed times are realized as the flow is in- 
creased to 200% of the anti-ice flow. 

Figure 4.3.2 is a sequence of photographs from ref. 1 showing 
the progressive shedding of a typical ice cap. Note that by the 
time the leading edge ice is shed, the runback of glycol along the 
wing surface has substantially removed the frost and ice particles 
behind the active portion of the leading edge. Operating in a deice 
mode may not be feasible if the aircraft can not tolerate the 
possible increase in drag and decrease in lift associated with the 
accreted ice prior to shedding. 
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Figure 4.3.2. - Deice Time Sequence for an Ice Accretion Formed at V = 96 Knots, LWC = 2.40 q/m } 
Angle of Attack = 7.8°, T = 23° F, and Icing Time = 10 Minutes. ’ 



4.4 Fluid Type Comparison 


Since two fluids are commercially available for use with the 
fluid ice protection system, it was felt desirable to test both 
fluids during part of the research to determine any differences in 
required flow rates for anti-icing. Both fluids, AL5 and TKS80, 
were tested on Wing Model A with the stainless steel mesh panel. 
Thirty four icing conditions, at 5 and 25 degrees F, were tested 
using both fluids. The fluid flow required for anti-icing using the 
AL5 fluid, SFFj. was compared with the fluid flow required when 
using the TKS80 fluid, SFT^, as presented in Table III. 

The results indicate that at the warm temperature, 25 degrees 
F, there is little difference between the two fluids over the 19 
runs compared at 96 and 145 knots. However, at the cold temperature 
of 5 degrees F, approximately 20 percent higher flow rate is re- 
quired when using the TKS80 fluid compared to the AL5 fluid. This 
occurs at both 96 and 145 knots. This difference cannot be directly 
explained by observing the freezing temperature profiles of the two 
fluids shown in Figure 2.8. Both fluids appear to have nearly iden- 
tical freezing temperature characteristics. 
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Table III: Comparison of Glycol Based Freezing Depressant Solutions 


V T 

Knots deg. F 


Runs Compared 
(runs) 


(SFF 2 - SFF^lOO 
(SFFj)(no. of runs) 


96 & 145 25 19 

96 & 145 5 15 


5.0 Anti-Ice Flow Rate Prediction Techniques 

5 . 1 ADS -4 Method 


Prior to 1983, the only published method of predicting anti-ice 
flow rates was that found in "The Engineering Summary of Airframe 

Icing Technical Data (ADS-4)" (ref. 6). This empirical method is 
based on the calculation of a glycol flow rate from the entire 
porous panel sufficient to lower the overall collection of water to 
a temperature equal to what is referred to as the datum temperature. 
This method does not take into account the distribution of the local 
impingement efficiency. The datum temperature has been used to 
represent the temperature of an unheated surface in icing. This 
method was not developed specifically for wing sections with porous 
leading edge panels. In fact, ref. 6 compares the predicted flow 
rates with experimental flow rates from wind tunnel tests using a 
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spray nozzle to distribute glycol on a r a dome (ref. 7). 

As mentioned, this empirical method is based an the calculation 
of the overall collection efficiency, E^. This collection is 
defined as the ratio of water collected to the maximum possible 
water catch. This is computed by the equation = (Y u -Y^)/h where 
Y^ is the starting ordinate of the upper water droplet tangent 
trajectory to the airfoil, Y^ the starting ordinate of the lower 
tangent trajectory, and h is the projected height of the airfoil. 

In reference 6 E^, Y^, and Y^ are graphically correlated with a 
dimensionless inertia parameter, Ko, for a variety of airfoil sec- 
tions. However, only a small portion of the airfoils in use today 
are included. Thus, a matching procedure is necessary to find an 
airfoil section that has been correlated that closely matches the 
airfoil of interest. It is often necessary to make extrapolations 
and assumptions when determining the values of E^, Y^, and Y y 
throughout the procedure listed in ref. 6. These approximations 
will in turn affect the accuracy of the predicted anti-ice flow 
rate . 

Since the publication date of ref. 6, computer codes have 
become available for computing limiting tangent trajectories, 
overall collection ef ficiencies , . and the local impingement efficien- 
cy distribution (Fig. 5.1). A program developed by Ohio State 
University under a grant from NASA LeRC (ref. 8) was used to cal- 
culate the values of E^, Y^, Y u> and h/c. This modified prediction 




method, as outlined below, will be referred to as the ADS-4 method. 
Appendix D contains a listing of the computer program used to cal** 
culate the ADS-4 predicted flow rates presented in this report. 

Step 1: Calculate the overall collection efficiency factor, E , as 

a function of 

(1) airfoil shape 

(2) airspeed 

(3) air density 

(4) drop diameter 

(5) angle of attack 
using a 2-D water droplet trajectory code. 

Step 2: Calculate the rate of water impingement on the airfoil from 

the equation 

(V) (LWC ) (E ) (h/ c ) 

M = .0031 (1) 

w (S - S. ) 

u 1 

3 2 

with LWC in g/m , V in knots, and in g/min cm . 

Step 3: Determine the glycol mass fraction, G, required to produce 

a solution with a freezing temperature equal to the datum tempera- 
ture. The datum temperature, as defined in ref. 6, represents the 
temperature of an unheated surface in icing, described as the "wet 
air boundary layer" temperature. This temperature is a function of 


airspeed, ambient temperature, and altitude. Reference 6 contains 
several plots at various altitudes shoving the relationship between 
the datum temperature with the ambient temperature and velocity. 


Otep 4: Calculate the fluid flow required to achieve the glycol 

mass fraction, 0, given the water catch rate, M y , by the equation 


SFFe 


(G)(M) 

_ w 
(X - G) 


( 2 ) 


where X is the initial glycol plus proponal mass fraction of the 
solution as it is pumped through the porous panel. 
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5.2 An alytical Method 


An analytical method of predicting anti-ice flow rates was 
first presented in reference 9. This method was developed after the 
original IRT tests in 1980 (ref. 1) of the fluid system using the 
stainless steel mesh panel on Wing Model B. A basic assumption of 
this method is that the minimum glycol fluid flow rate required for 
anti-icing results in a freezing temperature of the glycol water 
mixture equal to the average temperature between ambient and total 
temperature at the location of the maximum impingement efficiency. 
This average temperature is an arbitrary, empirical assumption ap- 
plicable when testing at velocities less than 175 knots. Further 
research is needed to determine if this assumption is valid at 
velocities above 175 knots. The local impingement efficiency (rate 
of water catch) varies in a chordwise direction, with the maximum 
efficiency located close to the stagnation location. This assump- 
tion was verified experimentally by noting that as the glycol flow 
rate is reduced from a flow rate in excess of that required for 
anti-icing, ice particles first begin to form at the stagnation 
location. However, at an angle of attack with a stagnation location 
close to the edge of the porous region this assumption is no longer 
valid. 

Any reasonably accurate 2-D water droplet trajectory program 
can be used to calculate the local maximum impingement efficiency, 
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8 However, the calculation of is not as straight forward 

max ' max 

as the calculation of E used in the ADS-4 method. The value of 

ID 

5 is determined from the maximum slope of a curve fit to the 
max 

final water drop trajectory positions, S and Y, on the airfoil as 
shown in Figures C.l - C.3 of Appendix C. These Figures present the 
trajectory positions and the local impingement efficiency distribu- 
tion for three different icing conditions with a different number of 
impacting trajectories at each condition. The range of for 

different numbers of trajectories varies from one icing condition to 
the next. 

The analytical method is presented here in a step-by-step 
procedure. Appendix D contains a listing of a computer program used 
to predict the analytical flow rates presented in this report. 

Step 1: Calculate the maximum local collection efficiency factor, 

6 as a function of; 
max 

(1) airfoil shape 

(2) airspeed 

(3) air density 

(4) drop diameter 

(5) angle of attack 

using any reasonably accurate 2-D water droplet trajectory computer 
code. 
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Step 2: Calculate the water catch rate, M# , by the equatio, 

M w = -0031*3^ (LWC)(V) 

With LWC in g/. J , V in knots, and M, i„ g/min cn 2 


Step 3= Determine the g!yco, mass fraction, G, required to 

a solution with a freezing temperature equal to the average 

the ambient and the stagnation temperatures. Normal!., this 

temperature is within 2 degrees F of the datum temperature 
the ADS-4 method. 


4: Calculate the fluid flow required to acnieve th. 

Pius proposal mass fraction, G, given a water catch rate, H 
equation ' 


SFFa = 


(G)(M ) 

w_ 

(X - G) 


"here X is the initial glycol plus proposal mass 
glycol solution that is supplied to the panel. 


fraction o 


( 3 ) 

produce 
between 
average 
used in 

glycol 
by the 

(4) 

f the 
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5.3 Comparison with Experimental Da ta 


Appendix B contains tables of the anti-ice averaged experimen- 
tal flow rates and the corresponding ADS-4 and analytical predicted 
flow rates. Included in the tables are the airfoil projected 
height, tangent trajectory limits, overall collection efficiency, 
and the maximum impingement efficiency. Predicted flow rates were 
computed for most, but not all the icing conditions. Also, only 
conditions at angles of attack between 1 and 6 degrees were com- 
puted. Table B.l presents the flow rates with the stainless steel 
panel on Wing Model A. Table B.2 presents the flow rates with the 
stainless steel panel on Wing Model B. Table B.3 presents the flow 
rates with the titanium panel on Wing Model A. Figure 5.2 presents 
the comparison between the experimental and ADS-4 predicted flow 
rates with (a) stainless steel panel on Wing Model A, (b) stainless 
steel panel on Wing Model B, and (c) titanium panel on Wing Model A. 
Figure 5.3 presents the same comparison between the experimental and 
analytical predicted flow rates. 
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Figure 5.2. - Comparison of Averaged Experimental and ADS-4 Predicted Anti-Ice 
Fluid Flow Rates. 
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(b) Wing Model B, Stainless Steel Mesh Panel 
Figure 5.2. - (Continued.) 
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(c) Wing Model A, Titanium Panel 

Figure 5.2. - (Continued.) 
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(a) Wing Model A, Stainless Steel Mesh Panel 

Figure 5.3. - Comparison of Averaged Experimental and Anal- 
ytically Predicted Anti-Ice Fluid Flow Rates. 
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b) Wing Model B, Stainless Steel Mesh Panel 


Figure 5.3. - (Continued.) 
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(c) Wing Model A, Titanium Panel 

Figure 5.3. - (Continued.) 
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To analyze the comparison statistically, a general purpose 
statistical computing program (Mini tab) was used. Table IV presents 
the results of the analysis. This analysis was done using all the 
flow rates except the two high flow rate conditions shown on Figures 
5.2 and 5.3 b and c, shown by the dashed symbols. These conditions, 
at V - 96 knots, T = 5 degrees F, LWC =2.4 g/m 3 , and d = 20 
n * crons * wet* tested only one time, and are significantly higher 
than the Continuous Maximum envelope conditions. It is felt that 
the large influence these points could have on the statistical 
analysis is not warranted by the uncertainty associated with these 
points . 

The correlation coefficient is an indication of the association 
between the experimental and predicted flow rates, with a value of 
1-0 representing a perfect correlation. The correlation coef- 
ficients are quite high, with the lowest coefficient of .852 . The 
highest correlation coefficient, .947, occurs with the stainless 
steel panel on the Wing Model A when comparing the analytical 
predicted to the experimental flow rates. Comparing the correlation 
coefficients for each wing-panel configuration, the analytical 
method is better correlated with the experimental results than the 
ADS-4 method. 

Also included in Table IV are the linear term coefficients 
determined through a least squares regression, while forcing the 
regression line th ough the origin (0,0). The predicted flow rate 
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was treated as the independent variable, and the experimental flow 
rate as the dependent variable as shown in equation 5. 

Coefficient*(Predicted Flow) = Experimental Flow (5) 

Thus, a coefficient of greater than 1.0 means the method predicts 
flow rates less than the actual flow rates. Opposing trends, as il- 
lustrated in Figures 5.2 and 5.3, are quickly noticed. The ADS-4 
method consistently predicts values less than the experimental flow 
rates, with coefficients ranging from 1.22 to 1.77. The analytical 
method consistently predicts values more conservative (higher) than 
the experimental flow rate, with coefficients ranging from .78 to 
.98, a much narrower range. The fact that the analytical method ap- 
pears to be more consistent from one configuration to the next is 
surprising since the calculation of S (jjax is less certain than the 
calculation of Em, used in the ADS-4 method. The fact that the 
analytical method tends to over predict the actual flow rates up to 
25 percent could be advantageous from a design and safety stand- 
point. The analytical method also has the important advantage of 
being more consistent from one wing-panel configuration to the next. 
It is felt that the analytical method is superior to the ads-4 
method for calculating anti- ice flow rates for a fluid ice protec- 
tion system. 
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Table IV: Coaparison of Experimental and Predicted Anti-Ice 

Fluid Flow Rates. 


ADS- 4 Method 



Linear Term 

Correlation 


Coefficient 

Coefficient 

Wing Model A, 



Stainless Steel Panel 

1.16 

.935 

Wing Model B, 



Stainless Steel Panel 

1.77 

.867 

Wing Model A, 



Titanium Panel 

1.22 

.880 


Analytical Method 



linear Term 

Correlation 


Coefficient 

Coefficient 

Wing Model A, 



Stainless Steel Panel 

.78 

.947 

Wing Model B, 



Stainless Steel Panel 

.98 

.931 

Wing Model A, 



Titanium Panel 

.86 

.934 


109 


5.4 Sensitivity Analysis 


The following examination of the possible areas of uncertainty 
in the data is necessary to appreciate the accuracy of the ex- 
perimental and predicted flow rates. 

(1) Accuracy of setting and measuring the glycol 

fluid flow rate, ± .001 g/miu cm 

(2) Variability of liquid water content in tunnel 
icing cloud, approximately r 20% 

(3) Variation of air temperature in tunnel test 
section, ± 1 degrees F 

(4) Variability in computing the overall collection 
efficiency E^, approximately ± 1% 

(5) Variability in computing the local 

impingement efficiency 8 , approximately t 10% 

(6) Repeatability of judging anti-ice flow rates, 

as high as 100%, but generally less than t 25 %. 

6.0 General Coaments 


Several comments can be made concerning the operation of the 
fluid ice protecti • and how the results can be applied to 
other situations. 
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1. Injecting a small amount of unfiltered glycol fluid through the 
titanium panel created a partial clogging of the porous region. 
This limited the amount of anti-ice data obtained at the end of the 
test plan, and did not allow reliable deice data to be acquired. A 
similar amount of unfiltered fluid did not affect the stainless 
steel panel. 

2. Often a flow rate of 25% lower than the anti-ice flow rate will 
still provide adequate protection while not allowing a permanent 
accretion of ice on the leading edge. 

3. Because of the subjectiveness in determining the actual anti-ice 
flow rates, the actual flow rate values should not be over em- 
phasized. Even the calculated averaged flow rates should be used 
cautiously since not all icing conditions were repeated enough to 
obtain a flow rate with a high degree of repeatability. 

A. A disadvantage of the icing conditions (LWC and d) available for 
testing in the IRT is that they do not directly correspond with the 
FAR Continuous Maximum and Intermittent Maximum Envelopes used for 
certification of ice protection systems. The glycol flow rates 
presented in this report, depending on the icing condition, may be 
several orders of magnitude higher than the flow rate for an icing 
condition within the FAR Part 25 Envelopes . 
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5. The flow rates presented in this report are typical of large 
leading edge radii general aviation airfoils only. The actual flow 
rates presented in this report cannot be generalized to other wing 
sections . 


7.0 Conclusions 


As a result of the tests presented in this report, the following 
conclusions have been made: 

1. Anti-ice flow rates are unaffected by angle of attack as long as 
the stagnation location is within the center 50 to 75 percent of the 
active porous region of the leading edge panel. 

2. As expected, the anti-ice flow rate increases with an increase 
in liquid water content and drop diameter. 

3. A new laser drilled titanium panel tested for the first time in 
the NASA LeRC IRT proved successful. No significant difference in 
required anti-ice flow rates between it and the stainless steel 
panel were noticed. 

4. The composite panel tested suffered from excessive porosity, 
causing the flow rates for anti-icing to be much higher than those 
required for the stainless steel panel. 

5. Deice shed times of three minutes or less are possible using a 
glycol fluid flow rate equal to the anti-ice flow rate. 

6. At 25 degrees F there is little difference between the TKS80 and 
AL5 glycol-water solutions. However, at 5 degrees F a higher flow 
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rate by approximately 20 percent is required when using the TKS80 
fluid compared with that required when using the AL5 fluid. 

7. An empirical method, ADS-4, predicts anti-ice flow rates less 
than the actual flow rates. This underprediction varied from one 
wing and panel configuration to the next. 

8. An analytical method developed under a NASA grant predicts flow 
rates up to 25 percent above the actual flow rates. This method was 

more consistent in its prediction for three different wing and panel 
configurations. 

8 • 0 Recommendations 


Additional icing wind tunnel tests ere planned using a titanium 
and a composite panel, ainilar to those discussed in this report, on 
a new natural laminar flow airfoil. Several recommendations for 

further evaluation of the fluid ice protection system are as fol- 
lows . 

1. Determine anti-ice flow rates at lower liquid water contents. 
These icing conditions are now possible due to recent modifications 
to the IRT spray nozzle system, and will more closely match the FAR 
Part 25 Continuous Maximum conditions. 
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2. Determine the best leading edge temperature to be used in the 
analytical prediction technique, especially at velocities over 200 
knots . 

3. Test the fluid ice protection system on an aircraft in natural 
icing. This may be accomplished using the NASA LeRC Twin-Otter 
Icing Aircraft or with a Cessna 206. The anti-ice, deice, and pump 
characteristics should be compared to those obtained during the wind 
tunnel tests. The anti-icing flow rates should be used to verify 
the prediction techniques to natural icing conditions . 

4. Further evaluate the deice characteristics through additional 
icing tunnel tests. Examine any differences in shed times depending 
if the system had been turned on temporarily prior to the ice accre- 
tion, with a minimal flow of glycol during the ice accretion, and 
with the system turned off prior to and during the ice accretion. 

5. Conduct wind tunnel tests with airfoils and airspeeds typical of 
high performace aircraft. 

6. Develop instrumentation to more objectively determine the 
minimum required anti-ice flow rates. 
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7. Examine the feasibility of a panel with chord wise varying 
porosity to direct a higher concentration of glycol to the stagna- 
tion region where the requirement for the glycol is the greatest. 

8. Compare the stagnation point locations determined experimentally 
to those predicted analytically. 
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Appendix A 




* 

l 


S n—ii ry of Anti - Ice Wind Tunnel Results 


Table A. 

1: Summary of Anti- 

-Ice Flow 

Rates on 

Wing Model 


A, 

Stainless Steel Mesn Panel, al5 

' Fluid. 

V 

T 

Alpha 

LWC 

d 


Knots 

°F 

degrees 

g/m J 

microns 

g/min cm* 

96 

25 

-0.5 

1.16 

11 

.0145 

96 

25 

2.1 

1.16 

11 

.0145 

96 

25 

2.1 

1.16 

11 

.0106 

96 

25 

5.9 

1.16 

11 

.0094 

96 

25 

7.7 

1.16 

11 

.0149 

96 

25 

2.1 

1.16 

13 

.0110 

96 

25 

2.1 

1.16 

13 

.0152 

96 

25 

5.9 

1.16 

13 

.0164 

96 

25 

5.9 

1.16 

13 

.0106 

96 

25 

5.9 

1.16 

13 

.0086 

96 

25 

-2.0 

1.16 

15 

.0489 

96 

25 

-0.5 

1.16 

15 

.0196 

96 

25 

2.1 

1.16 

15 

.0133 

96 

25 

2.1 

1.16 

15 

.0133 

96 

25 

2.1 

1.16 

15 

.0184 

96 

25 

5.9 

1.16 

15 

.0135 

96 

25 

5.9 

1.16 

15 

.0115 

96 

25 

7.7 

1.16 

15 

.0239 

96 

25 

9.5 

1.16 

15 

.0346 

96 

25 

2.1 

1.50 

15 

.0156 

96 

25 

2.1 

1.50 

15 

.0172 

96 

25 

2.1 

1.50 

15 

.0153 

96 

25 

5.9 

1.50 

15 

.0188 

96 

25 

5.9 

1.50 

15 

.0160 

96 

25 

5.9 

1.50 

15 

.0125 

96 

25 

9.5 

1.16 

11 

.0669 

96 

25 

2.1 

1.80 

15 

.0196 

96 

25 

2.1 

1.80 

15 

.0172 

96 

25 

5.9 

1.80 

15 

.0156 

96 

25 

5.9 

1.80 

15 

.0215 

96 

25 

-0.5 

1.55 

20 

.0297 

96 

25 

-0.5 

1.55 

20 

.0239 

96 

25 

0.5 

1.55 

20 

.0203 

96 

25 

2.1 

1.55 

20 

.0160 

96 

25 

2.1 

1.55 

20 

.0258 

96 

25 

5.9 

1.55 

20 

.0180 

96 

25 

5.9 

1.55 

20 

.0215 

96 

25 

5.9 

1.55 

20 

.0219 

96 

25 

5.9 

1.55 

20 

.0160 

96 

25 

7.7 

1.55 

20 

.0786 

96 

25 

7.7 

1.55 

20 

.0368 

96 

25 

9.5 

1.55 

20 

.0223 

96 

25 

2.1 

2.40 

20 

.0246 

96 

25 

5.9 

2.40 

20 

.0305 
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f 



V 

T 

Alpha 

LWC 

d 

SFF’x 

Knots 

°F 

degrees 

g / m 3 

microns 

g/min 

96 

25 

5.9 

2.40 

20 

.0321 

96 

25 

5.9 

2.40 

20 

.0215 

96 

5 

- 0.5 

1.16 

11 

.0332 

96 

5 

- 0.5 

1.16 

11 

.0309 

96 

5 

2.1 

1.16 

11 

.0372 

96 

5 

2.1 

1.16 

11 

.0282 

96 

5 

5.9 

1.16 

11 

.0458 

96 

5 

5.9 

1.16 

11 

.0270 

96 

5 

5.9 

1.16 

11 

.0254 

96 

5 

7.7 

1.16 

11 

.0317 

96 

5 

7.7 

1.16 

11 

.0258 

96 

5 

9.5 

1.16 

11 

.0645 

96 

5 

5.9 

1.16 

13 

.0407 

96 

5 

5.9 

1.16 

13 

.0387 

96 

5 

5.9 

1.16 

13 

.0282 

96 

5 

- 0.5 

1.16 

15 

.0571 

96 

5 

- 0.5 

1.16 

15 

.0415 

96 

5 

- 0.5 

1.16 

15 

.0630 

96 

5 

- 0.5 

1.16 

15 

.0461 

96 

5 

2. 1 

• 1.16 

15 

.0618 

96 

5 

2.1 

1.16 

15 

.0407 

96 

5 * 

2.1 

1.16 

15 

.0477 

96 

5 

2.1 

1.16 

15 

.0403 

96 

5 

5.9 

1.16 

15 

.0587 

96 

5 

5.9 

1.16 

15 

.0461 

96 

5 

5.9 

1.16 

15 

.0352 

96 

5 

5.9 

1.16 

15 

.0403 

96 

5 

7.7 

1.16 

15 

.0540 

96 

5 

7.7 

1.16 

15 

.0280 

96 

5 

7.7 

1.16 

15 

.0540 

96 

5 

7.7 

1.16 

15 

.0501 

96 

5 

2.1 

1.50 

15 

.0512 

96 

5 

5.9 

1.50 

15 

.0610 

96 

5 

5.9 

1.50 

15 

.0567 

96 

5 

5.9 

1.50 

15 

.0391 

96 

5 

2.1 

1.80 

15 

.0559 

96 

5 

5.9 

1.80 

15 

.0774 

96 

5 

2.1 

1.55 

20 

.0763 

96 

5 

2.1 

1.55 

20 

.0786 

96 

5 

5.9 

1.55 

20 

.0696 

96 

5 

7.7 

1.55 

20 

.0829 

96 

5 

2.1 

2.40 

20 

.1189 

96 

5 

5.9 

2.40 

20 

.0962 

96 

5 

5.9 

2.40 

20 

.0481 

96 

5 

5.9 

2.40 

20 

.0919 

96 

5 

5.9 

2.40 

20 

.0915 

96 

-10 

5.9 

1.16 

11 

.0598 
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V 

T 

Alpha 

LWC 

d 

SFF X 

Knots 

°F 

degrees 

g / m 3 

microns 

g/min cm 2 

96 

-10 

5.9 

1.16 

15 

.0888 

96 

-10 

5.9 

1.50 

15 

.1146 

145 

25 

- 0.5 

0.80 

11 

.0168 

145 

25 

2.1 

0.80 

11 

.0133 

145 

25 

2.1 

0.80 

11 

.0149 

145 

25 

5.9 

0.80 

11 

.0117 

145 

25 

7.7 

0.80 

11 

.0375 

145 

25 

7.7 

0.80 

11 

.0395 

145 

25 

9.5 

0.80 

11 

.0946 

145 

25 

2.1 

0.80 

15 

.0242 

. 145 

25 

2.1 

0.80 

15 

.0282 

145 

25 

5.9 

0.80 

15 

.0203 

145 

25 

5.9 

0.80 

15 

.0336 

145 

25 

- 0.5 

1.16 

15 

.0774 

145 

25 

- 0.5 

1.16 

15 

.0418 

145 

25 

0.5 

1.16 

15 

.0336 

145 

25 

2.1 

1.16 

15 

.0262 

145 

25 

2.1 

1.16 

15 

.0289 

145 

25 

2.1 

1.16 

15 

.0270 

145 

25 

5.9 

1.16 

15 

.0340 

145 

25 

5.9 

1.16 

15 

.0286 

145 

25 

7.7 

1.16 

15 

.0650 

145 

25 

7.7 

1.16 

15 

.0571 

145 

25 

7.7 

1.16 

15 

.0594 

145 

25 

2.1 

1.16 

19 

.0442 

145 

25 

5.9 

1.16 

19 

.0430 

145 

25 

5.9 

1.16 

19 

.0383 

145 

25 

2.1 

1.55 

20 

.0493 

145 

25 

2.1 

1.55 

20 

.0489 

145 

25 

5.9 

1.55 

20 

.0540 

145 

25 

7.7 

1.55 

20 

.0391 

145 

25 

7.7 

1.55 

20 

.0544 

145 

25 

2.1 

1.70 

20 

.0415 

145 

25 

5.9 

1.70 

20 

.0505 

145 

5 

- 0.5 

0.80 

11 

.0798 

145 

5 

2.1 

0.80 

11 

.0704 

145 

5 

5.9 

0.80 

11 

.0508 

145 

5 

5.9 

0.80 

11 

.0469 

145 

5 

5.9 

0.80 

11 

.0579 

145 

5 

2.1 

0.80 

15 

.0743 

145 

5 

5.9 

0.80 

15 

.0681 

145 

5 

5.9 

0.80 

15 

.0763 

145 

5 

- 0.5 

1.16 

15 

.1134 

145 

5 

2.1 

1.16 

15 

.1048 

145 

5 

2.1 

1.16 

15 

.1134 

145 

5 

5.9 

1.16 

15 

.0915 

145 

5 

5.9 

1.16 

15 

.0860 
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V 

T 

Alpha 

Knots 

°F 

degrees 

145 

5 

7.7 

145 

5 

7.7 

145 

5 

2.1 

145 

5 

5.9 

145 

5 

2.1 

145 

5 

2.1 

145 

5 

5.9 

145 

5 

5.9 

145 

5 

7.7 

145 

5 

2.1 

145 

5 

5.9 

145 

-10 

5.9 

145 

-10 

5.9 

145 

-10 

5.9 


LWC 

g/m 3 

d 

microns 

SFF X 

g/min cm 

1.16 

15 

.0892 

1.16 

15 

.1033 

1.16 

19 

.1044 

1.16 

19 

.1384 

1.55 

20 

.1867 

1.55 

20 

.1862 

1.55 

20 

.1560 

1.55 

20 

.1474 

1.55 

20 

.1760 

1.70 

20 

.1392 

1.70 

20 

.1224 

0.80 

11 

.0857 

0.80 

15 

.1115 

1.16 

15 

.1295 
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Table A. 2: Summary of Anti-Ice Flow Rates on Wing Model 

A, Stainless Steel Mesh Panel, TKS8Q Fluid. 


V 

T 

Alpha 

Knots 

°F 

degrees 

96 

25 

2.1 

96 

25 

5.9 

96 

25 

2. 1 

96 

25 

5.9 

96 

25 

2.1 

96 

25 

5.9 

96 

25 

5.9 

96 

25 

5.9 

96 

25 

2.1 

96 

25 

5.9 

96 

25 

5.9 

96 

25 

5.9 

96 

25 

2.1 

96 

25 

2.1 

96 

25 

5.9 

96 

25 

5.9 

96 

5 

2.1 

96 

5 

5.9 

96 

5 

2.1 

96 

5 

5.9 

96 

5 

2.1 

96 

5 

5.9 

96 

5 

5.9 

96 

5 

2.1 

96 

5 

5.9 

96 

5 

5.9 

96 

-10 

5.9 

96 

-10 

5.9 

96 

-10 

5.9 

96 

-10 

5.9 

145 

25 

2.1 

145 

25 

2.1 

145 

25 

5.9 

145 

25 

5.9 

145 

25 

5.9 

145 

25 

2.1 

145 

25 

5.9 

145 

25 

5.9 

145 

25 

5.9 

145 

25 

2.1 

145 

25 

2.1 

145 

25 

5.9 

145 

25 

5.9 

145 

25 

2 1 


LWC 

d 

SFFy 

g/m 3 

microns 

g/min cm 

1.16 

a 

.0117 

1.16 

ii 

.0090 

1.16 

15 

.0176 

1.16 

15 

.0145 

1.50 

15 

.0211 

1.50 

15 

.0368 

1.50 

15 

.0297 

1.50 

15 

.0196 

1.80 

15 

.0227 

1.30 

15 

.0278 

1.80 

15 

.0168 

1.80 

15 

.0336 

1.55 

20 

.0258 

1.55 

20 

.0278 

1.55 

20 

.0203 

1.55 

20 

.0254 

1.16 

11 

.0375 

1.16 

11 

.0372 

1.16 

15 

.0583 

1.16 

15 

.0438 

1.50 

15 

.0489 

1.50 

15 

.0559 

1.50 

15 

.0415 

1.80 

15 

.0626 

1.80 

15 

.0532 

1.80 

15 

.0579 

1.16 

15 

.1205 

1.16 

15 

.0845 

1.50 

15 

.1228 

1.50 

15 

.0849 

0.80 

11 

.0211 

0.80 

11 

.0203 

0.80 

11 

.0203 

0.80 

11 

.0235 

0.80 

11 

.0235 

0.80 

15 

.0254 

0.80 

15 

.0223 

0.80 

15 

.0282 

0.80 

15 

.0282 

1.16 

15 

.0344 

1.16 

15 

.0360 

1.16 

15 

.0375 

1.16 

15 

.0293 

1.55 

20 

.0508 


122 


V 

T 

Alpha 

LWC 

d 

SFFy 

Knots 

°F 

degrees 

g/m 3 

microns 

g/min 

145 

25 

5.9 

1.55 

20 

.0684 

145 

25 

2.1 

1.70 

20 

.0516 

145 

25 

5.9 

1.70 

20 

.0360 

145 

5 

2.1 

0.80 

11 

.0763 

145 

5 

2.1 

0.80 

11 

.0571 

145 

5 

5.9 

0.80 

11 

.0618 

145 

5 

5.9 

0.80 

15 

.0626 

145 

5 

5.9 

0.80 

15 

.0724 

145 

5 

2.1 

1.16 

15 

.1119 

145 

5 

2.1 

1.16 

15 

.1087 

145 

5 

5.9 

1.16 

15 

.0829 

145 

5 

5.9 

1.16 

15 

.1138 

145 

5 

2.1 

1.55 

20 

.1572 

145 

5 

5.9 

1.55 

20 

.1435 

145 

-10 

2.1 

0.80 

11 

.1095 

145 

-10 

5.9 

0.80 

11 

.0845 

145 

-10 

5.9 

0.80 

11 

.0927 

145 

-10 

5.9 

1.16 

15 

. 1416 
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Table A. 3: Summary of Anti-Ice Flow Rates on Wing Model 

8, Stainless Steel Mesh Panel, Observer 1. 


V 

T 

Alpna 

LWC 

d 

SFF! 

Knots 

°F 

degrees 

g/m 3 

microns 

g/min cm 2 

96 

25 

7.8 

1.16 

11 

.0068 

96 

25 

4.5 

1.16 

11 

.0055 

96 

25 

1.2 

1.16 

U 

.0086 

96 

25 

12.0 

1.16 

15 

.0126 

96 

25 

7.8 

1.16 

15 

.0102 

96 

25 

4.5 

1.16 

15 

.0138 

96 

15 

1.2 

1.16 

15 

.0108 

96 

25 

4.5 

1.55 

15 

.0138 

96 

25 

4.5 

1.80 

15 

.0194 

96 

25 

4.5 

1.55 

20 

.0200 

96 

25 

7.8 

2.40 

20 

.0305 

96 

25 

4.5 

2.40 

20 

.0272 

96 

25 

1.2 

2.40 

20 

.0317 

96 

9 

4.5 

2.17 

17 

.0690 

96 

9 

4.5 

2.40 

20 

.0805 

96 

5 

12-0 

1.16 

11 

.0487 

96 

5 

7.8 

1.16 

11 

.0253 

96 

5 

4.5 

1.16 

11 

.0349 

96 

5 

1.2 

1.16 

11 

.0302 

96 

5 

4.5 

2.40 

20 

.0869 

96 

5 

1.2 

2.40 

20 

.0940 

145 

25 

1.2 

1.16 

15 

.0388 

175 

25 

7.8 

0.80 

11 

.0346 

175 

25 

4.5 

0.80 

11 

.0355 

175 

25 

1.2 

0.80 

11 

.0376 

175 

25 

7.8 

1.16 

15 

.0635 
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Table A. A: 


Summary of Anti-Ice Flow Rates on Wing Model 
0, Stainless Steel Mesh Panel, Observer 2. 


V 

T 

Alpha 

LWC 

d 

sff 2 

Knots 

°F 

degrees 

g/m 3 

microns 

g/min cm 2 

96 

25 

7.8 

1.16 

11 

.0120 

96 

25 

4.5 

1.16 

11 

.0117 

96 

25 

12.0 

1.16 

15 

.0346 

96 

25 

7.8 

1.16 

15 

.0269 

96 

25 

4.5 

1.16 

15 

.0219 

96 

25 

4.5 

1.16 

15 

.0135 

96 

25 

4.5 

1.55 

15 

.0182 

96 

25 

4.5 

1.80 

15 

.0185 

96 

25 

4.5 

1.55 

20 

.0426 

96 

25 

4.5 

1.55 

20 

.0200 

96 

25 

1.2 

1.55 

20 

.0416 

96 

25 

1.2 

1.55 

20 

.0247 

96 

25 

7.8 

2.40 

20 

.0292 

96 

25 

4.5 

2.40 

20 

.0336 

96 

25 

1.2 

2.40 

20 

.0394 

96 

5 

7.8 

1.16 

11 

.0200 

96 

5 

7.8 

1.16 

11 

.0450 

96 

5 

4.5 

1.16 

11 

.0370 

96 

5 

1.2 

1.16 

11 

.0379 

96 

5 

7.8 

1.16 

15 

.0484 

96 

5 

7.8 

1.16 

15 

.0475 

96 

5 

4.5 

1.16 

15 

.0379 

96 

5 

4.5 

1.16 

15 

.0577 

96 

5 

4.5 

1.16 

15 

.0465 

96 

5 

1.2 

1.16 

15 

.0426 

96 

5 

1.2 

1.16 

15 

.0699 

96 

5 

4.5 

1.55 

15 

.0435 

96 

5 

4.5 

1.80 

15 

.0715 

96 

5 

4.5 

1.55 

20 

.1097 

145 

25 

1.2 

1.16 

15 

.0542 

145 

25 

1.2 

1.16 

15 

.0265 

145 

25 

1.2 

1.55 

20 

.0629 

175 

25 

7.8 

0.80 

11 

.0209 

175 

25 

4.5 

0.80 

11 

.0256 

175 

25 

1.2 

0.80 

11 

.0327 

175 

25 

4.5 

0.80 

15 

.0443 

175 

25 

7.8 

1.16 

15 

.0404 

175 

25 

4.5 

1.16 

15 

.0349 

175 

25 

4.5 

1.16 

15 

.0567 

175 

25 

1.2 

1.16 

15 

.0479 

175 

25 

4.5 

1.16 

19 

.0675 

175 

25 

4.5 

1.16 

19 

.0706 

175 

25 

1.2 

1.55 

20 

.0641 

175 

25 

1.2 

1.55 

20 

.0817 
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Table A. 5: Summary of Anti-Ice Flow Rates on wing Model 

A, Titanium Panel, Observer 1. 


V 

T 

?• ipha 

nots 

°F 

degrees 

96 

25 

-2.0 

96 

25 

-.5 

96 

25 

2.1 

96 

25 

5.9 

96 

25 

7.7 

96 

25 

2.1 

96 

25 

5.9 

96 

25 

-.5 

96 

25 

2.1 

96 

25 

5.9 

96 

25 

,7.7 

96 

5 

-2.0 

96 

5 

-.5 

96 

5 

2.1 

96 

5 

5-9 

96 

5 

7.7 

96 

5 

2.1 

96 

5 

5.9 

96 

5 

-.5 

96 

5 

2.1 

96 

5 

2.1 

145 

25 

2.1 

145 

25 

5.9 

145 

25 

2.1 

145 

25 

2.1 

145 

25 

5.9 

145 

25 

2.1 

145 

25 

2.1 

145 

25 

5.9 

145 

25 

2.1 

145 

25 

2.1 

145 

25 

2.1 

145 

25 

2.1 

145 

5 

2.1 

145 

5 

5.9 


LWC 

d 

SFFi 

g/m 3 

macrons 

g/min c 

1.16 

15 

.0204 

1.16 

15 

.0194 

1.16 

15 

.0146 

1.16 

15 

.0094 

1.16 

15 

.0090 

1.50 

15 

.0166 

1.50 

15 

.0152 

1.55 

20 

.0329 

1.55 

20 

.0281 

1.55 

20 

.0211 

1.55 

20 

.0211 

1.16 

15 

.0745 

1. 16 

15 

.0759 

1.16 

15 

.0648 

1.16 

15 

.0551 

1.16 

15 

.0492 

1.50 

15 

.0672 

1.50 

15 

.0575 

1.55 

20 

.0735 

1.55 

20 

.0710 

1.80 

15 

.0679 

0.80 

15 

.0281 

0.80 

15 

.0225 

1.16 

15 

.0388 

1.16 

15 

.0340 

1.16 

15 

.0291 

1.16 

19 

.0440 

1.16 

19 

.0402 

1.16 

19 

.0423 

1.55 

20 

.0627 

1.55 

20 

.0541 

1.70 

20 

.0683 

1.70 

20 

.0582 

0.80 

15 

.0665 

0.80 

15 

.0617 

• 


Table A. 6: Summary of Anti-Ice Flow Rates on Wing Model 

A, Titanium Panel, Observer 2. 


V 

T 

Alpha 

LWC 

d 

sff 2 

Knots 

°F 

degrees 

g/m 3 

microns 

g/min cm 2 

96 

25 

-2.0 

1.16 

15 

.0374 

96 

25 

-2.0 

1.16 

15 

.0166 

96 

25 

-.5 

1.16 

15 

.0159 

96 

25 

2.1 

1.16 

15 

.0152 

96 

25 

5.9 

1.16 

15 

.0101 

96 

25 

7.7 

1.16 

15 

.0076 

96 

25 

2.1 

1.50 

15 

.0159 

96 

25 

5.9 

1.50 

15 

.0104 

96 

25 

-.5 

1.55 

20 

.0280 

96 

25 

2.1 

1.55 

20 

.0198 

96 

25 

2.1 

1.55 

20 

.0218 

96 

25 

5.9 

1.55 

20 

.0184 

96 

25 

7.7 

1.55 

20 

.0180 

96 

25 

2.1 

1.80 

15 

.0121 

96 

25 

2.1 

2.40 

20 

.0256 

96 

25 

5.9 

2.40 

20 

.0222 

96 

5 

-.5 

1.16 

15 

.0520 

96 

5 

2.1 

1.16 

15 

.0490 

96 

5 

2.1 

1.16 

15 

.0534 

96 

5 

5.9 

1.16 

15 

.0336 

96 

5 

5.9 

1.16 

15 

.0460 

96 

5 

7.7 

1.16 

15 

.0492 

96 

5 

7.7 

1.16 

15 

.0500 

96 

5 

2.1 

1.50 

15 

.0530 

96 

5 

5.9 

1.50 

15 

.0515 

96 

5 

2.1 

1.55 

20 

.0730 

96 

5 

5.9 

1.55 

20 

.0710 

96 

5 

2.1 

1.80 

15 

.0621 

96 

5 

2.1 

2.40 

20 

.0841 

96 

5 

5.9 

2.40 

20 

.0770 

145 

25 

2.1 

0.80 

15 

.0381 

145 

25 

5.9 

0.80 

15 

.0277 

145 

25 

5.9 

1.16 

19 

.0444 
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Table B.2: Averaged Experimental and Predicted Flow Rates on 

Winq Model B with Stainless Steel Mesh Panel. 
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Table B.3: Averaged Experimental and Predicted Flow Rates on 

Wing Model A with Titanium Panel. 
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Water Droplet Impingement Efficiency Distributions 



Figure C.l. - Local Impingement Efficiency Distributions at V = 96 Knots 
T = 25° F, Angle of Attack = 5.9°, and d = 13 microns. 






Figure C.l. - (Continued. 






Impingement Efficiency Distributions 
0 F, Angle of Attack = 2.1°, and d = 



(Continued. ) 



Figure C.2. - (Continued. ) 





(a) With 4 Trajectories Impacting Airfoil 

Figure C.3. - Local Impingement Efficiency Distributions at V = 96 Knots 
T = 25° F, Angle of Attack = 2.1°, and d = 20 microns. 



Figure C.3. - (Continued.) 



Appendix D 


^ gg uter — °S raM Listing of Prediction Methods 

100 REM ADS* 4 Prediction of Minimum Required Anti-Ice Flow Rates 
110 REM for a Fluid Ice Protection System 
120 REM Enter the Following Values 
130 PRINT "Wing Chord Length ?"; 

140 INPUT C 

150 PRINT "liquid Water Content, LWC (g/ra*3) ?"; 

160 INPUT LWC 

170 PRINT "True Velocity, Vt (Knots) ?"; 

180 INPUT VKTS 

190 PRINT "Datum Temperature, Tok (deg. F) ?" ; 

200 INPUT TDAT 

210 PRINT "Overall Water Collection Efficiency, Em ?"; 

220 INPUT Em 

230 PRINT "Which Glycol Solution (1 = AL5 or 2 =TKS80) ?"; 

240 INPUT FLUID 

250 PRINT "Chordwise Impingement Extent, Su - SI ?"; 

260 INPUT SUSL 

270 PRINT "Airfoil Projected Height, h/c ?"; 

280 INPUT HC 
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290 REM Computation Section 

300 Ren Conpute Overall Water Catch Rate 

310 MW = . 379*VTKS*LWC*C*HC*EM/ (SUSL - C) 

320 IF FLUID = 1 THEN X = 87. 8: G = 36.698-.69*TDAT-.014' ,f TDAT’ v TDAT 

330 IF FLUID * 2 THEN X = 8l.2:G = 36.565-.67*TDAT-.013*TDAT*TDAT 

340 REM Compute Predicted Anti-Ice Fluid Flow 

350 SFFE = ((MW*G)/(X - G))*. 00856 

360 PRINT "Predicted Anti-Ice Fluid Flow ";SFFE 

370 STOP 

380 END 
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Appendix D.2 




* 


Listing of Computer Program in BASIC of 
Analytical Prediction Method 

100 REM Analytical Prediction of Minimum Required Anti-Ice 
110 REM Flow Rates for a Fluid Ice Protection System 
120 REM Enter the following values 
130 PRINT ’’Wing Chord Length, C (ft) ?"; 

140 INPUT C 

150 PRINT "Liquid Water Content, LWC (g/a"3) ?"; 

160 INPUT LWC 

170 PRINT "True Velocity, Vt (Knots) ?"; 

180 INPUT VKTS 

190 PRINT "Total Temperature, T (deg. F) ?"; 

200 INPUT TT0T 

210 PRINT "Max. Local Water Collection Efficiency, Bmax ?" 
220 INPUT BETAMAX 

230 PRINT "Which Glycol Solution (1 = AL5 or 2= TKS80) ?" ; 

240 INPUT FLUID 

250 REM Computation Section 

260 REM Convert velocity to cm/min 

270 VCM = VKTS*3088 . 84 

280 REM Compute local water catch rate 


290 WATER = VCM*LVC*BETAMAX/100~3 

300 REM Compute Average Temperature 

310 TSTAT = TTOT-((VCM/1828.82)'2./12020) 

320 TAVG = (TSTAT+TT0T)/2 

330 IF FLUID = 1 THEN X = 87. 8: G = 36.698-.69*TAVG-.014*TAVG*TAVG 
340 IF FLUID = 2 THEN X = 81. 2: G = 36.565-.67*TAVG-.013*TAVG*TAVG 
350 REM Compute Predicted Anti-Ice Fluid Flow 
360 SFFA = (WATER*G)(X-G) 

370 PRINT "Predicted Fluid Flow = ";SFFA 
380 STOP 
390 END 
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